Stacking Sequence Dependence of Graphene Layers on SiC(OOOl) - Experimental and 

Theoretical Investigation 
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Different stacking sequences of graphene are investigated a using combination of experimental 
and theoretical methods. The high-resolution transmission electron microscopy (HRTEM) of the 
stacking sequence of several layers of graphene, formed on the C-terminated 4H-SiC(0001) surface, 
was used to determine the stacking sequence and the interlayer distances. These data prove that the 
three metastable configurations exist: ABAB, AAAA, ABCA. In accordance to these findings, those 
three cases were considered theoretically, using Density Functional Theory calculations comparing 
graphene sheets, freestanding and positioned on the SiC(OOOl) substrate. The total energies were 
calculated, the most stable structure was identified and the electronic band structure was obtained. 
The four graphene layer electron band structure depends crucially on the stacking: for the ABAB 
and ABCA stacking, the bands, close to the K point, are characterized by the hyperbolic dispersion 
relation while the AA stacking the dispersion in this region is linear, similar to that of a single 
graphene layer. It was also shown that the linear dispersion relation is preserved in the presence of 
the SiC substrate, and also for different distances between adjacent carbon layers. 

PACS numbers: 61.50.Ah, 81.10.Aj 
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Graphene is one of the most extensively investigated 
perspective materials at the moment. Its unique elec- 
tronic properties, promising very important applications 
in electronics, have attracted interests of many scientific 
groups. The massless fermion dispersion, observed in 
the vicinity of the K point in the energy range close 
to the Fermi surface, and the related, experimentally 
determined the ^/B dependence of the Landau levels, 
was reported from investigations of a single, freestand- 
ing graphene layer [H-f^. In order to implement graphen 
in the electronics industry, a method for the synthesis 
of reproducible, good quality graphene layers, stabilized 
on a solid substrate, has yet to be developed. The most 
natural one is the deposition of the epitaxial graphene 
on the SiC substrate, either on the Si- or C-terminated 
principal faces of the SiC crystal which is currently con- 
sidered as one of the most promising method to attain 
this goal A recent work Q shows that the method 

is potentially able to furnish good quality devices. 

Graphene multilayer structures have been studied both 
experimentally and theoretically. From these studies, it is 
recognized that the multilayer graphene electronic prop- 
erties are strongly dependent on the stacking sequence. 
Both the Bernal (AB) and rhombohedral (ABC) stack- 
ing sequences are described theoretically [lO, EH and ex- 
perimentally [13, [3] . It was also shown that for rhom- 
bohedral(ABC) stacking the band gap can be induced 



by an external electric field [lj|. The description of the 



AA structure was often overlooked because it is energet- 
ically unfavourable |l5|. However, very recently, Nori- 
matsu and Kusunoki [13| have observed the existence of 
the AA stacking sequence around the step on the Si- 



face. We have undertaken systematic investigations of 
the four layer graphene synthesized on the C-terminated 
4H-SiC(0001) surface, such as the one which has been 
investigated recently |16j. By an optical inspection, the 
three different regions were identified, which are then 
subjected to the HRTEM scrutiny. The results are pre- 
sented in Fig. [TJ The three different stacking sequences 
are identified, showing the existence of the AA stack- 
ing sequence as well as the Bernal (AB) and rhombohe- 
dral (ABC) sequences, independent of the presence of the 
steps. Thus these three structures could be synthesized 
on the C-terminated 4H-SiC(0001) surface proving that 
these structures are metastable, and could be obtained 
by the appropriate selection of the thermodynamic pa- 
rameters of SiC annealing Q. These structures could 
persist for macroscopically long times. 

The HRTEM investigations were also used to deter- 
mine the distances between the graphene sheets. The 
example of such investigation, shown in Fig. [l] presents 
the determination of the carbon interlayer distances. The 
high precision measurements prove that the carbon bot- 
tommost layer is separated by a relatively large distance 
of 3.0 A from the SiC surface, which is close to the inter- 
layer spacing in graphite, the situation which is a stan- 
dard for graphene on the C-terminated face [3]. There- 
fore these graphene layers are not strongly bound to the 
underlying SiC surface, in contrast to the layers grown 
on the opposite, Si-terminated surface The following 
carbon sheet is separated by a distance of 3.7A, while 
the remaining one by 3.4A, i.e. essentially equal to the 
carbon sheet separation in graphite. These experimen- 
tally determined graphene interlayer distances were used 
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FIG. 1: (Colour online) (left) A cross-sectional HRTEM image of graphene layers , grown on the C-terminated face of 4H-SiC. 
The first graphene layer is separated by 3.0A from the SiC surface, next - by 3.7A from the previous one, and the subsequent 
layers are separated equally by 3.4A . (right) Experimental HRTEM filtered images showing different stacking sequences of 
graphene layers. The positions of the carbon atoms columns are marked by circles. 



in the DFT ab initio calculations described below. In the 
present work, we employ ab initio the density functional 
theory the VASP [17-20j code to investigate graphene- 
SiC interface. We have used the projector augmented 
wave (PAW) approach 21] in its variant available in the 
VASP package [20] . The local spin density approximation 
(LSDA) was applied for the exchange-correlation func- 
tional. The plane wave cutoff energy was set to 500 eV. 
The Monkhorst-Pack k-point mesh was set to 7 x 7 x 1. 
The 4H-SiC(0001) supercell was constructed using 8 bi- 
layers of Si-C. Four graphene layers were located at the 
top of the SiC(OOOl) surface, at the separation deter- 
mined experimentally (see Fig. [1]). Since the graphene 
interlayer distance results from the van der Waals in- 
teraction, which, as a rule, could not be obtained from 
DFT calculations properly, the combined experimental- 
theoretical approach is the only possible way to obtain 
the properties of multilayer graphene with high precision. 
In the present calculations, the slab replicas are separated 
by the space of about 17 A. An elastic adjustment was 
performed at the interface due to the lattice mismatch 
between SiC and graphite such that the two top SiC lay- 
ers and the graphene layer was relaxed in the plane. The 
conjugate gradient algorithm was used in the relaxation 
of the atomic positions. The model, a/S x •\/3i?30° — SiC 
unit cell with a fitted graphene layer (GL) (22| - [23 l. was 
used. The direct DFT calculation indicate that for free- 
standing graphene layers, Bernal (ABAB) stacking se- 
quence is the most stable, having the total energy equal 
to -80.824 which is lower than the energy of both 
the rhombohedral (ABC A) -80.836 eV and (AAAA) - 
80.782 eV stackings, which is consistent with the earlier 
works [31. The electron energy bands are analyzed us- 



ing a comparison between the results obtained for four 
graphene layers, freestanding and deposited on the 411- 
SiC(OOOl) surface as presented in Fig. [2] The three 
stacking configurations, AAAA, ABAB and ABCA ex- 
hibit the miscellaneous electronic structure due to the 
different symmetry. For the case of the AAAA config- 
uration, the Dirac-type dispersion relation is observed. 
This is in agreement with a simple, tight binding argu- 
ment presented recently by Gonzalez et al [2^. As ex- 
pected, the band structure is different for both the Bernal 
(ABAB) and rhombohedral (ABCA) structures. In the 
case of the Bernal (ABAB) structure, in accordance to 
the earlier results, both the conduction and valence bands 
are hyperbolic [l^. The zero bandgap structure has the 
Fermi energy located at the common maximum of the 
valence and the minimum of the conduction bands. A 
slightly different picture is obtained for the freestanding 
rhombohedral (ABCA) sequence. The hyperbolic depen- 
dence, similar to that observed for the Bernal stacking 
is preserved. The two bands constitute a zero bandgap 
structure at the K point. Other bands are collagsed at 
the K point in the so-called "wizard-hat" shape [12]. Due 
to the lower crystallographic symmetry of the rhombo- 
hedral structure, they have the energy minimum shifted 
away from the K point (see Fig. [2]). 

The influence of the 4H-SiC substrate is different in 
these three cases. In the case of the (AAAA) stack- 
ing, the presence of the SiC substrate amounts to mere 
addition of the additional energy bands, of their ener- 
gies located in the valence band energies of graphene. 
In the case of Bernal stacking, the influence of the sub- 
strate leads to a slight modiflcation of the band structure, 
close to the K point. Thus the epitaxial graphene on 
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FIG. 2; (Colour online) Left - the spatial distribution of electron density, right - band structures of four graphene layer system 
in different stacking sequences: AAAA (top), ABAB (middle), ABCA (bottom); the left column - freestanding graphene, the 
right column - graphene deposited on the SiC(OOOl) substrate. 



the SiC(OOOl) surface is characterized by the essentially 
identical band structure for both (AAAA) and Bernal 
(ABAB), in respect to the freestanding graphene in these 
configurations, that indicates a relatively weak coupling 
of the graphene with the SiC surface as suggested by the 
relatively large distance measured by the HRTEM (see 
Fig. [1]). In contrast to that, in the case of rhombohe- 
dral stacking, the structure is most affected, leading to 
the opening of the gap, and the creation of the double 
minimum in the conduction band. The second conduc- 
tion band is concave and it is not on the Fermi level (see. 
Fig. [2|). The first valence band which is curved up at the 
K point, is associated with the SiC substrate influence. 

In order to explain why the symmetry of the systems 
so strongly affects the electronic structure, we have calcu- 
lated electron densities between the first two freestand- 
ing graphene layers in AAAA, ABAB, ABCA stacking 
sequence Fig. [5] (left). In order to compare the electron 
density distribution in these three cases, we have used a 
consistent scale. For the AAAA stacking sequence, the 
areas of the increased electron density (red pattern) are 
arranged in the honeycomb pattern. This electron dis- 
tribution pattern is similar to the isolated one graphene 
layer. Indeed, such a "highway" for electrons results in 
linear dispersion bands. The situation is quite different 
for graphene layers arranged into ABAB and ABCA se- 
quences. Electron distributions for these two cases are 
similar and they form isolated islands of the increased 
electron density (red pattern). This different electron 
configuration gives hyperbolic HOMO and LUMO bands 



dispersion (see Fig [2] (left) ) . 

In order to determine the properties of the graphene 
AA stacking sequence, the bi-, three- and four-layer free- 
standing graphene structure was simulated using DFT 
calculations. The resulting band structure is shown in 
Fig 131 From these results, it follows that the band struc- 
ture is qualitatively similar for a various number of layers. 
For all cases, the dispersion relation is linear, with the 
shift arising from the overlap of the adjacent carbon lay- 
ers. For a pair of carbon layers, the overlap of pz orbitals 
is changed by a change of the interlayer distance (see Fig. 
[3]- top), resulting in a different mutual shift of the two 
bands. As shown in Fig. [3] (bottom), for an odd number 
of graphene layers, one of the Dirac cone is exactly at 
the K point and is crossed by the Fermi level as in the 
isolated graphene layer. For an even number of layers, 
the intersection of the bands is located symmetrically in 
the vicinity of the K point. 

In this letter, using the HRTEM images, we have 
demonstrated that the epitaxial graphene, grown on 
the SiC substrate, can exist in three metastable con- 
figurations, i.e. the Bernal (ABAB), the rhombohedral 
(ABCA) and the AAAA stacking sequences. It was 
shown that the AAAA stacking, by virtue of its sym- 
metry, has pronouncedly different electronic properties, 
being linear in the vicinity of the K point. The other 
two configurations have their properties different adopt- 
ing the hyperbolic dispersion relation in the vicinity of 
the K point. The two structures are affected by the pres- 
ence of the SiC substrate. In contrast to that, the AAAA 
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FIG. 3: (Colour online) The band structure for two, three and 
four freestanding graphene layers (bottom) . Top - the depen- 
dence of AE shifting of the two bands of the AA stacking on 
the distance between carbon layers. 

stacking graphene preserves its linear dispersion relations 
in the presence of the SiC substrate. It is also demon- 
strated that the change of the distance between carbon 
sheets amounts to a mere shift of the band, still preserv- 
ing its linear character. Thus, the obtained results open 
the route to mechanically stable, fast electronic devices, 
using the Dirac dispersion relation typical for the AAAA 
stacking. 

This profound difference in electronic properties be- 
tween AAAA and other ABCA and ABAB stacking was 
confirmed by the direct plots of electron densities. It was 
observed that the charge pattern for the AAAA stack- 
ing is similar to an isolated graphene layer. It was also 
shown that for ABAB and ABCA structures charge accu- 
mulation regions are arranged into isolated islands. It is 
expected that the honeycomb pattern is more favourable 
for electron traffic than the insular one which confirms 
the advantage of the AAAA graphene as a material for 
future fast electronic devices. 
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